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The neoclassical polarization current, generated when a magnetic island propagates 
O 
CN through a tokamak plasma, is believed to influence the initial stage of the neoclassical 

C3 tearing mode evolution. Understanding the strength of its contribution in the relevant 
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plasma collision frequency regimes for future tokamaks such as ITER is crucial for 
the successful control and/or avoidance of the neoclassical tearing mode. A nonlinear 
drift kinetic theory is employed to determine the full collision frequency dependence 
Ch of the neoclassical polarization current in the small island limit, comparable to the 

j^ trapped ion orbit width. Focusing on the region away from the island separatrix 

^ (where a layer with a complex mix of physics processes exists) , we evaluate for the first 

_ O time the variation of the neoclassical ion polarization current in the transition regime 

>-> between the analytically tractable collisionless and collisional limits. In addition, 

Qh the island propagation frequency-dependence of the neoclassical polarization current 
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and its contribution to the island evolution is revealed. For a range of propagation 
frequencies, we find that the neoclassical polarization current is a maximum in the 
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Q^^ intermediate collision frequency regime analyzed here - a new and unexpected result 
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I. INTRODUCTION 

A tokamak plasma is subject to a number of magnetohydrodynamic (MHD) instabilities 
that can limit the performance of the toroidal fusion device. In tearing mode instabilities, 
filamentation of the plasma current density along equilibrium magnetic field lines forms a 
chain of magnetic islands at a rational surface, and leads to a corrugation of flux surfaces in 
its vicinity. Within the island, the enhanced radial transport of particles and heat reduces 
the radial pressure gradient. This results in a reduction of the core plasma pressure, which 
degrades the plasma confinement. According to the theory of single-fiuid resistive MHD, the 
evolution of a magnetic island is characterized by the rate of change in its width, which is 
proportional to the parameter A -: a measure of the free energy stored in the plasma current 
density for the magnetic reconnection to occur. For A' > 0, the island is predicted to grow. 
The neoclassical theory of tearing modes incorporates the effects of toroidal geometry in 
the layer surrounding the rational surface. One contribution to the current that influences 
the island evolution is the bootstrap current, which is proportional to the radial pressure 
gradient. Because of the pressure gradient flattening, the bootstrap current is suppressed 
inside the island region. This perturbation in the bootstrap current enhances the original 
fllamentation of the plasma current and hence drives the island growtlP^. The neoclassical 
tearing modes (NTMs) are characterized by this enhanced drive for the island growth, whose 
strength diminishes with increasing island width, w. The result is that w typically saturates 
at a substantial fraction of the tokamak minor radius, r. In toroidal geometry, the curva- 
ture of the magnetic fleld lines provides a stabilizing contribution to the island evolution 
equation"^ ^. This curvature effect is also found to be proportional to the pressure gradient 
and diminishes with the island width. Thus it can be thought of as an effect that reduces the 
bootstrap drive (except for sufficiently small islands^, though the effect is typically small 
for the large aspect ratio tokamaks we consider in this paper. 

There is substantial experimental evidence for the existence of a threshold mechanismP^^^, 
whereby a sufficiently small "seed" island (typically of O(lcm)) does not grow to a large 
saturated island, but rather shrinks away. However, there is no concrete theoretical frame- 
work to explain the observed level of the threshold width for the island width and provide 
quantitative predictions. One possible mechanism is the effect of finite cross-field transport 
of particles and heat in the vicinity of the island separatrixpHSl. A consequence of the finite 



radial transport is that the pressure gradient is not completely flattened across the island 
width, and then the bootstrap drive for the island width is reduced. Another possible source 
of the threshold mechanism is the polarization current, which is induced when the island is in 
relative motion with the bulk of the plasma. A number of works have considered the flnite 
ion Larmor radius (FLR) effect on tearing mode evolution in sheared slab geometry^^Hia^ 
which can give rise to the polarization current when the island width is comparable to the 
ion Larmor radius, pn. However, whether this FLR effect can stabilize the island depends 
on the plasma conditions and, in particular, the island propagation frequency. In toroidal 
geometry, the effect of finite trapped ion orbit width, or the ion "banana width", p;,j, can also 
generate a polarization current. Because pu ^ pui this neoclassical contribution dominates 
the polarization current except perhaps in the vicinity of the island separatrix, where the 
FLR effect is also likely to be important. The origin of the neoclassical polarization current 
can be understood as follows. As the magnetic island propagates through the plasma, ions 
and electrons respond differently; the ion response is dominated by the E x B drift, whereas 
the parallel transport dominates the electron response. However, because of the difference 
in ion and electron banana widths {pu ^ Pbe), the orbit-averaged E x B drifts of trapped 
ions and electrons differ, resulting in a net current perpendicular to the magnetic field lines 
(see Fig. [I]). This is the neoclassical polarization current, which in turn generates a parallel 
current perturbation to ensure that V.J = 0, where J is the current density. It is this paral- 
lel current perturbation that contributes to the island evolution. Past analytical works^^"^ 
show that this neoclassical polarization current contribution depends strongly on the plasma 
collision frequency regime; it is 0(e^/^) smaller in the collisionless limit {vu <^ ew), compared 
to the collisional limit {vu ^ too). Here, vu is the ion-ion collision frequency, e = r/R is 
the inverse aspect ratio (r and R are the minor and major radii respectively) and u is the 
island propagation frequency in the E x B rest frame. In this paper, we aim to develop 
the nonlinear drift kinetic theory of Refs. 1201 and |2I] further to determine the full collision 
frequency dependence of the neoclassical polarization current. 

As discussed above, the neoclassical polarization current is seeded by the trapped particles 
(through the finite ion banana width effect). Collisions transfer the current perturbation 
carried by the trapped particles to the passing particles at the rate vu/e. In the high collision 
frequency limit vu/t ^ uj (but collisionality, v^ <1), the collisional momentum transfer takes 
place sufficiently quickly that the time variation of the trapped particle response is resolved, 
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FIG. 1. A cross section of a magnetic island indicating the directions of the electric field pertur- 
bation (E), equilibrium magnetic field (Bq), island propagation [uj) and the polarization current 



and the polarization current is communicated to the passing particles and amplified in the 
process. On the other hand, in the low collision frequency limit {va <^ eu) the collisional 
transfer does not resolve the time variation and the polarization current remains at the seed 
level, which is 0(e^/^) smaller than that in the collisional limit. 

Taking into account the contributions discussed so far, the island evolution equation can 
be written in the form of the "modified Rutherford equation" : 
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Here, r^j is the resistive time scale, A;,^ denotes the bootstrap drive for the island growth 
and Aggj denotes the contribution from the curvature effect. As discussed in detail below, 
we separate the contribution from the neoclassical polarization current into two parts: Apoi 
denotes the contribution from the neoclassical polarization current away from the island 
separatrix (which we shall call the "external polarization current") and Aiayer denotes the 
contribution from the narrow layer that surrounds the island separatrix (the "layer polariza- 
tion current"). Our interest in this paper is focused on the term Apoi. We characterize the 
strength of the neoclassical polarization current by introducing the dimensionless parameter 
g{i'ii, e, u) and write 



\oi = -5'(«^M,e,w) 
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where Lq = q/{dq/dr), Lp = p/{dp/dr) and /3g is the plasma poloidal /3. Analytic results 



from Ref . |20] in the coUisionless limit and Refs. [19] and [21] in the colhsional hmit provide 

J lMe^^^f{iyii/euj)uj{u - uj^pi)/ujli vu/eu < 1, 
g{vii,e,uj) = < (3) 

y uj{uj - uj^pi - kT]iUj^i)/u^. Vii/euj > 1, 

where /(uu/eu) = 1 + ^^Vajtuj and w*pj = w*j(l + r\i). Here, w^,,, = {mTi/qqin){dn/dx) is 
the ion diamagnetic frequency, k = —1.17 in the banana regime^, ?7i = Ln/Lxi is the ratio of 
the density to temperature gradient length scales, L„ = n/{dn/dr) and Lti = Ti/{dTi/dr). 



7 is a weak logarithmic function of vu/euj and the dependence on '^^valeuj appears from 
the finite colhsional effects in the narrow layer around the trapped/passing boundary in 
pitch angle space^^. Other works have considered the effect of collision frequency on the 
polarization current, including: Ref. [231 which employs the linear MHD theory; and Refs. 
[Ml and [251 in which the collision frequency dependence is calculated by numerical simulation. 
In Refs. [2j|and[25l however, only the ion response is determined and quasineutrality is not 
imposed to derive the self-consistent electrostatic potential, $. As has been pointed out in 
Refs. [201 and [211 the polarization current is proportional to the second radial derivative of 
$, $". Hence it is likely to be important to work with a self-consistent form of $, which 
introduces an additional complication. 

There is a large contribution (formally a (5-function) to $" from the vicinity of the island 
separatrix, where simplified models predict a discontinuity in the electric field; this generates 
a skin current contribution to the ion polarization current. This provides an additional 
contribution to the island evolution equation denoted by lS.iayer in Eq. ([I]), which is expected 
to be comparable in magnitude and of opposite sign to the contribution from outside the 
separatrix layer^^'^, Apoi. However, the contribution from this separatrix layer requires an 
extended physics model to treat it, including the FLR effeclP', non-linearities in the parallel 
electric field and non-perturbative cross-field diffusion^^. Consequently its contribution to 
the island evolution scales differently with plasma parameters to that outside the layer. 
It is therefore appropriate to separate the calculation into two regions, considering the 
contribution to the neoclassical ion polarization current from each region separately. In this 
paper, we address only the contribution from outside the separatrix layer, which we called 
the "external" polarization current [see Eq. ([I|]. This enables us to focus on the physics of 
the transition between the collision frequency regimes without the complicated physics of 
the separatrix layer. We leave the more complicated layer polarization current calculation 



for future research. 

To summarize, the objective of this paper is to extend the drift kinetic theory developed 
in Refs. [SU] and [5T] to determine the full collision frequency dependence of the "external" 
neoclassical ion polarization current (i.e. that away from the separatrix layer) and calcu- 
late its contribution to the island evolution. In particular, we consider the dependence of 
(yf(i/jj, e, oj) on the collision frequency in the previously unexplored intermediate regime, where 
Vii/euj ~ 1. In addition, we consider the influence of the island propagation frequency u on 
the contribution of the neoclassical polarization current to the island evolution. We find that 
whether Apoi provides a stabilizing contribution (in this paper, this corresponds to (7 > 0) 
depends crucially on the relative size of u compared to the ion diamagnetic frequency, u^^i. 

The paper is organized as follows: in Section|Tl], we summarize the calculations of particle 
responses to the perturbed magnetic geometry, which is discussed in full in Ref. [20] At 
leading order in pu/w, we introduce the free function associated with integration along 
magnetic field lines, hi, which carries the leading order collision frequency dependence due to 
the finite ion banana width effect. Constraint equations at higher order in pu/w determine 



its full form. In Section III we discuss the method of calculating the collision frequency 



dependence giyu^e^uj) from the ion response to the magnetic island. Our new results for 



g{i'ii,e,uj) across the collision frequency range are presented in Section IV, and conclusions 
are drawn in Section |Vl 

II. MAGNETIC GEOMETRY AND THE DRIFT KINETIC EQUATION 

A tokamak plasma with a large aspect ratio (e <^ 1) and a circular cross section is 
considered. We introduce a single helicity magnetic perturbation of the following form, 
assuming the constant-?/' approximation: 

^(e)=^cose, (4) 

where ^ is the helical angle: 

,^,„(«-i). (5) 

Here, g^ = Tn/n is the value of the safety factor at the rational surface, and m and n are the 
poloidal and toroidal mode numbers respectively. Then, in the toroidal coordinate system 
(x, 0, 0), where x is the poloidal flux, 9 is the poloidal angle and is the toroidal angle, the 



total magnetic field is given by: 

B = J(x)V0 + V0x V(x + ^), (6) 

where I{x) = RB^ and B^ is the toroidal component of B. It is convenient to define a 
perturbed flux function satisfying B.Vfi = 0: 

fi=^^^^^-cose, (7) 

where w^ = RBqw, Bg is the poloidal field and Xs is the value of x at the rational surface. 
The particle responses to the magnetic island perturbation are described by the drift 
kinetic equation: 
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where v^; is the E x B drift, v^ = —v\\h x V(f n/wcj) is the magnetic drift and Ucj = qjB/rrij 
is the gyrofrequency. $ is the electrostatic potential perturbation to be obtained from 
quasineutrality, and C is the model collision operator. In the coordinate system {x,0,O^ 
the parallel derivative operator is given by: 
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(9) 



and q'^ = {dq/dx)x=xs- -^^ discussed in the introduction, we restrict our analysis to the 
region outside the separatrix layer, which corresponds to fi > 1. The parallel derivative 
operator can then be annihilated by introducing the following averaging operators: 



and 



{■■■)o = ^f-d9, (12) 

(■■■)e = E ^ / •••^^' (13) 
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where Eqs. ( 12 ) and ( 13 ) are the 6'-averaging operators for the passing and trapped particles 
respectively, and 9^ is the bounce point for the trapped particles. In Eq. (13), a = ±1 is 
the sign of the parallel velocity: f y = avy/l — XB. Vy in Eq. (N) can be annihilated by the 
operator {{Rq...)^)^. 

The perturbed distribution function in Eq. ([s]) is expressed in terms of the adiabatic and 
non-adiabatic parts: 



fj 



q^ 



Fmj + gj, 



(14) 



where Fmj is the Maxwellian distribution. The non-adiabatic term, Qj, is solved for each 
particle species by expanding it in terms of two small parameters^Sl A = w/r and Sj = Pbj/ui, 
where pbj = ^^^'^Pej is the trapped particle banana width and poi = niiVthi/qiBe- Thus, 

^. = E'^r^"^r"^' (15) 

and we solve for the expansion terms gj^'""^ by considering the relevant order contributions 
to the drift kinetic equation ([s]). For the purpose of this paper, we are interested in the 
collision frequency dependence of the ion response. The full derivations of the ion and 
electron responses are given in Ref. EOl Here, we restrict our discussion to a description of 
the key parts of the calculation. 

To leading order (0((5fA°)), the electron and ion distribution functions are given by 
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where the bar above a quantity indicates that it is independent of 6*, 
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and the self-consistent electrostatic potential is 
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Here, h{VL) is a free function that is related to the electron density profile in the vicinity of 
the island separatrix and can be determined from the consideration of the radial particle 
transportpSl^. The 0{5l^'^) equation provides: 
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where hi{Q, $,) is a free function that arises as a consequence of integration along unperturbed 
field lines. It is shown in Ref. [201 that the leading order contribution to the collision 
frequency dependence of the polarization current conies from this term. The explicit form 
for hi is determined from a solubility constraint on the higher order equation. This constraint 
equation is obtained by averaging the 0{6iA) equation over the unperturbed field lines. The 
equation for the passing particles in the limit u ^ A;||t'|| (appropriate for thin islands) is: 
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while for the trapped particles: 
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where |f||| = v\\/l — \B\. In the coUisionless limit ua <^ ecu, the second term of Eq. (21) [or 



Eq. (22)] is negligible and the first term, which describes the response to the E x B flow, 



determines h^. In the opposite limit ua ^ ecu, the first term of Eq. (21) becomes negligible 



and the collisional effect alone determines hi. These analytic limits have been explored in 
Refs. E3 and ED In this paper, we consider the arbitrary collision frequency regime between 



the two analytically tractable limits and numerically solve the full constraint equations (21) 



and (22). In our calculation we use the following momentum-conserving model collision 
operator^ 



Ciiig-i) = 2vii{v) 
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for the ion-ion collisions. Here, z/jj(f) is the 90° deflection frequencySSl 
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where x = v/vthi, 4>{x) is the error function, G{x) is the Chandrasekhar function 

4>{x) — X(f)'{x) 



G{x) 



2x2 



(25) 



and 0' = d(f)/dx. In Eq. (23), the parallel flow in the momentum conservation term is 
defined as 

^¥ = — P — T^T / ^^^ ^ii{v)m9i, (26) 

where 



e "" /''^hiv^Vii{v) dv. 



Ion-electron collisions are small compared to the ion-ion collisions and are therefore ne- 



glected. In the next section, we discuss the method of solving Eqs.(21) and (22), and how 
the coefficient g{vii^e^u) is determined from the solution for hi. 



III. CALCULATION OF g{yii,e,u) 



As shown in Ref. [201 the ion response provides the dominant contribution to the piece of 
the parallel current perturbation, J\\ , which varies along magnetic field lines. This provides 
the contribution of the neoclassical ion polarization current to the island evolution. The 
leading order contribution to this parallel current perturbation comes from gi ' (which 
includes the collision frequency dependent part, hi): 
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Integrating Eq. (27) provides two parts for Jn: first is the parallel current perturbation 



arising from the neoclassical polarization current, which we identify as the part of J\\ that 
varies along the field lines and fiux surface- averages to zercP^. The second part, which is 
constant on a fiux surface, is the fiux surface average of the bootstrap current perturbation, 
which is determined from the perturbed ion and electron parallel fiows. It is incorporated 



in the function of Vt following the integration of Eq. (27). In this paper we focus on the first 



part: the contribution from the neoclassical polarization current, which is determined from 
Eq. (27) with the condition that this part of J\\ fiux surface-averages to zero. Once J\\ is 



determined, the contribution of the neoclassical polarization current to the island evolution 
can be calculated via 
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FIG. 2. Plots of the numerical results for hi as a function \B for different values of vu/eoj, at 
e = 0.1: (a) is at uu/eu) = 10^^; (b) at vu/euj = 10^^; (c) at vu/eu) = 10^^ and (d) at vu/euj = 10*^. 
The vertical dotted lines correspond to the trapped/passing boundary, A = Ac- The dashed curves 
in (b) ~ (d) correspond to the collisionless result, vajeuj = 10~^, to aid comparison. Note the 
steep gradients that need to be resolved in the vicinity of the trapped/passing boundary at low 
values of vu/euj. 

where the summation is over the region x> Xs and x < Xs- As discussed in Section |H the 
colhsion frequency dependence of the neoclassical polarization current and its contribution 
to the island evolution is described by the coefficient gi^ua^e^uj). 

A numerical code is developed to determine giyu^e^uj) from the solution for hi, using 
Eqs.(20),(21), (27), (28) and ([2]). A particularly careful treatment is required for analyzing 
the pitch angle space; as discussed in Ref. EU in the total absence of collisions (i.e. when the 



second term of Eq. (21 ) becomes zero), hi is discontinuous at the trapped/passing boundary. 
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Inclusion of collisional effects smooths out this discontinuity in a narrow boundary layer (see 
Fig. |2|. The pitch angle mesh needs to be closely packed to resolve this "dissipation layer" 
that surrounds the trapped/passing boundary, in which the collisional effect is important 
even in the low collision frequency limit, due to steep gradients in pitch angle. A further 
complication is the treatment of the flow term, uiu in the momentum conserving term of 



the model collision operator (see Eq. (23)). This provides an integro-differential equation, 
which we solve by an iterative numerical scheme. The results for gl^n, e, u) are presented in 
the next section. 

IV. RESULTS 

The numerical result for g{uii,e,uj) as a function of ua/ecu is shown in Fig. Isj As dis- 
cussed in Section [T], the convention in this paper is that g > corresponds to a stabilizing 
contribution to the island evolution. As expected, g(uii,e,u}) is 0{e^^'^) smaller in the coUi- 
sionless limit compared to the collisional limit, and agrees well with the analytic results in 
both limits. The transition from one limit to another takes place approximately between 
Uii/eco ~ 0.1 and vu/euj ~ 100, for e = 0.1. Compared to the previous linear MHD results^, 
the new results presented here show that 5'(z/jj, e, u) starts increasing from a somewhat lower 
value of Vii/euj] in the fitting made to the linear MHD theory in Ref. [23] g stays at the 
low collisionless value until vu/euj ~ 1. This difference is due to the ^Jva/tuo dependence of 
g{vii, e, oj) arising from the leading order collisional correction in the low collision frequency 
limit, as predicted by Ref. |2I1 This ^Jva/tuj dependence was omitted in the fitting sug- 
gested by Ref. |23l However, our results suggest that it is important even when ua/euj ^ 1, 
enhancing the collisionless neoclassical polarization current by a factor ~ 2 at vu/eu ~ 1, 
for example (compared to the collisionless value). In Fig. Ill we show that the leading order 



collisional correction to g in the collisionless limit is indeed 0{y^Uii/eijj)^ aside from the weak 
logarithmic dependence which offsets the gradient from the expected value of 1/2. This is 
consistent with analytic theory^^, which predicts such a correction arising from a narrow 
layer in pitch angle space around the trapped/passing boundary. 

Fig. [3] shows a case where g{h'ii,e,uj) is positive (and therefore the external polarization 
current provides a stabilizing contribution) across all of the collision frequency domain. It 
turns out, however, that the sign of g{i'ii, e, u) depends on the relative size of cj with respect 
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FIG. 3. Plot of giya^ e, w) as a function of collision frequency regime vn/euj for e = 0.1, r/, = 1.0 
and oj/uj^i = 2.5. The horizontal dotted lines represent the analytic values of g{i'ii, €,uj) in the two 
collision frequency limits. The plot shows that g starts deviating from its collisionless value from 
as early as uu/euj ~ 0.1, and this deviation is clear by uu/euj = 1. 
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FIG. 4. Plot of the gradient: d\og{IS.po\{i^ii) — Apoi(t'M = 0))/di^ii against uu/eu) for e = 0.1, 
r]i = 0.5 and oo/ujn = 2.5. The dashed line represents the gradient taken from the analytic result in 
Ref. [2n while the solid line shows the gradient taken from our new numerical result for g{i'ii, e,uj). 
If g{i'ii,e,u}) scaled as ^vajeiij, then the gradient is expected to be 1/2. The deviation from 1/2 
is due to the weak logarithmic dependence on va/eoj, predicted analyticall}^. 



to the ion diamagnetic frequency, a;*,. The analytic forms of g^va.e.u) in the collisionless 
and collisional limits [see Eq. ([s])] demonstrate that g is negative if oj/uj^i is positive and 
less than (1 + r/j) in the collisionless limit, or if it is less than 1 + (1 + k)rii in the collisional 
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FIG. 5. Plots of g{vii,e,oj) vs. vu/euj for e = 0.1 and rji = 1.0 for different values of w. The 
numbers on the right hand side of the graph are the values of a;/cj*j for each of the plots. 
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FIG. 6. Plots of (/(fjj, e, w) vs. ua/euj for e = 0.1 and r/j = 1.0 for smaller values of w, showing that 
g is minimum at uj/oj-^pi = 1/2 in the collisionless limit, consistent with the analytic form given in 
Eq. ([3]). The numbers on the right hand side of the graph are the values of ujIuj^i for each of the 
plots. 



limit. Fig. p^ shows how g varies with uj /u^i and collision frequency. For sufficiently small 
uj (< w^Kj), g is negative everywhere in the collision frequency domain, and it is positive 
everywhere for uj > aj^,j(l + T]i), as expected from Eq. ^. In Fig. |6| with e = 0.1 and 
rji = 1.0, we show that the minimum for g{vii,e,uj) in the collisionless limit is indeed at 
(jj/uj^i = 1, as expected from the dependence uj[uj — a;*j(l + rji)], with rji = 1.0. Again, we 
see that our numerical results are in agreement with the analytic results in the appropriate 
limits. 
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A particularly interesting feature that is evident from Fig. |5]is that the sign oi g changes in 
the intermediate collision frequency regime as uu/eu is increased, when Uint < ou ^ u;*i(l+r7j), 
where for the parameters used in Fig. [5} Umt — l-loo^i- Furthermore, there is a range of 
u where g has a maximum at intermediate collision frequencies between the two analytic 
limits. The implication of this is that whether the neoclassical polarization current can 
stabilize or amplify the magnetic island depends not only on plasma parameters, including 
the collision frequency regime, but also on the relative rotation between the island and 
the plasma. Conversely, whether or not the external polarization current heals or amplifies 
magnetic islands can depend sensitively on the collision frequency regime. In Fig. [7] we show 
that this behavior of g is robust for different values of e and rji. 

Finally, we consider the relationship between the coefficient k{e) and the critical value 
of oj in the coUisional limit Uc, for which the sign of g{vii,e,uj) reverses. According to Eq. 
PI, this critical uj is expected to be: ujc = uj*pi + k{e)7]iUJn. Using k{e) = — 1.17/c(e)P^ 
and explicitly calculating the passing particle fraction, /c, we find k = —0.67 and hence 
Uc = 1.33ci;=Kj, for e = 0.1 and rji = 1.0. This is consistent with our numerical result, where 
we see from Figs. |5] and |6] that Uc = I.SSw^kj. Lowering e to 0.02 we find k = —0.93 and have 
Uc = 1.03 with rji = 0.5, which is again in agreement with Uc obtained from the numerical 
solution of Eq. (21), shown in Fig. [TJ An important point to make is that in both cases, 
k{e) is substantially different from the infinite aspect ratio limit: ^(e = 0) = —1.17, even 
though e is very small. The significance of this is that the /c-dependence of k should be 
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properly taken into account in any quantitative analysis, even for a very small value of the 
inverse aspect ratio. 

V. CONCLUSION 

We have separated the contribution of the neoclassical polarization current to the island 
evolution into two parts: the "external" polarization current which exists outside the island, 
and a "layer" polarization current which exists in the complex boundary layer in the vicinity 
of the island separatrix. We have focused on the external contribution and determined 
the full collision frequency dependence of its contribution to the island evolution, using 
nonlinear drift kinetic theory. Our numerical results show that the collisional correction to 
this external contribution to the neoclassical polarization current is important even at very 
low collision frequency, Va/too ~ 0.1 (e = 0.1). Furthermore, we have found a rich structure 
in the contribution of the polarization current to the NTM evolution in the intermediate 
collision frequency regime which is not accessible by analytic theory. We find that whether 
it can provide a stabilizing contribution to the island evolution depends crucially on the 
size of oj/ujn as well as the collision frequency regime, vajtuj. Our results suggest that an 
element of an NTM avoidance or control scheme in future devices, such as ITER, may be 
through control of the plasma collision frequency as well as the rotation of magnetic islands. 
Future work will address the layer contribution to the polarization current. This opposes 
the external contribution that we discussed here, and therefore is important to include in 
order to make specific quantitative statements concerning the stabilizing influence of total 
neoclassical ion polarization current on NTMs. 
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